NPS ARCHIVE 
1966 
LEWIS, J. 




ANALYSIS AND SYNTHESIS OF A PURE FLUID 
STATE ATTITUDE CONTROL SYSTEM 

by 

Jack W# Lewis 

Thesis Supervisor: Herbert H. Richardson 

Date Submitted: May 20, 1966 






za©©@(jt>os[i©®‘ 

GENUINE PRESSBOARD BINDER 

cat. no. BF 2507 EMB 



ACCO OGDENSBURG, N.Y. 

CHICAGO, TORONTO, 

LONDON, MEXICO. D. F. 




ANALYSIS AND SYNTHESIS OF A FUSE FLUID 

* 

STATE ATTITUDE CONTROL SYSTEM 

V 

by 

Jaok W* Lewis 

// 



SUBMITTED IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE DEGREE OF 
NAVAL ENGINEER 
AND THE DEGREE OF 

MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

/ 

at the 



MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
May, 1966 



W' r f'tfcu >/£ 

(c? Co 



ANALYSIS AND SYNTHESIS OF A PURE FLUID 
STATE ATTITUDE CONTROL SYSTEM 



rl'y 

* > School 

u iy. Ciilfo^nia 



by 



Jaok W. Lewis 



i . 

Submitted to the Department of Naval Architecture and Karine 
Engineering on Kay 20, 1966, in partial fulfillment of the re- 
quirements for the Master of Scienoe degree in Meohenioal 
Engineering and the professional degree, Naval Engineer. 



ABSTRACT 

A pneumatio, pure fluid oontrol system, capable of main- 
taining the angular position of an inertia load in a gravita- 
tional space, using only one moving part, was designed and 
built. The control system consisted of a vortex-sink angular 
rate sensor, a pendulum controlled flapper-nozzle angular posi- 
tion sensor, a summing amplifier and bistable reaction jets. 

The system was designed to operate in a bang-ban^ manner in 
order to reduce errors in a minimum amount of tine. 

A previous author's work on the vortex-sink angular rate 
sensor is extended to include information and data on propor- 
tional amplification of the output signal and the frequency 
response of the sensor. Methods for synthesizing the components 
into a working system are presented along with a theoretical 
analysis of the system's performance. 

Results of experiments performed on the completed system 
showed that the reaction jets could be switohed in a nearly op- 
timal manner when -responding to errors introduced on one side 
of the system, but, due to the unsymmetrical gain of the rate 
sensor, the errors introduced on the opposite side resulted in 
a certain amount of after end-point chattering. 

t 

Good correlation between theory and experiment led to the 
conclusion that the application o^* standard oontrol system design 
techniques can be readily applied to the design of pure fluid oon- 
trol systems. 



Thesis Supervisor: Herbert H. Richardson 

Title: Associate Professor of Mechanical Engineering 



ii 



ACKNGV/LEDGMENTo 



I wish to express my sincere thanks to frofessor 
Richardson who first introduced me to the subjeot of pure 
fluid amplification and whose keen insight into the funda- 
mentals of dynamic systems and fluid control helped me over- 
come the many problems encountered in this work. I wish, 
also to mention that because of his lively ability as a 
professor, I first beoarae interested in the field of control 
systems and that without the warm professor-student relation- 
ship experienced with him, 1 might never have been able to 
accomplish this work. 

My thanks also to Messrs. James Sitomer, Stanley 
Goodwin and Paul Hinckley of the Instrumentation Laboratory 
who assisted me in gathering the data on the rate sensor. 

Special thanks are due Mr. William Ribich who sug- 
gested the design of the position sensor used in this work 
and who is responsible for the competent photography con- 
tained herein. 

I would like to thank Messrs. Adam Bell, Pa\il Blaiklook 
Chester Nachtigal, Richard Sidell and David Wormley of the 
Design ana Controls Division of the Mechanical Engineering 
Department for their friendliness and advice. 

Mrs. Martha Buckley is responsible for the competent 
typing of this thesis. •< 

This thesis was supported in part by the U.S. Air Foroe 
under oontraot AF 33(615)-2210 and sponsored by the Division 
of Sponsored Researoh of M.I.T. 



TABLE OF CONTENTS 



Abstract 11 

Acknowledgments . Ill 

Table of Contents lv 

Table of Figures • • « vl 

Chapter 

I. INTRODUCTION 1 

1.1 General Introduction to Pure Fluid State 

Technology ... .....1 

1.2 Purpose and Scope of Thesis 1 

1.3 Discussion of the Proposed System to Attain 

the Design Objectives 2 

n. VORTEX ANGULAR RATE SENSOR ANALYSIS AND DESIGN .... 5 

2.1 General Theory of Operation . . . 5 

2.2 Rate Sensor Design and Design Procedure. .’.... 20 

2.3 Signal Amplification 22 

2.4 Rate Sensor Dynamics 32 

2.5 Testing Procedure and Results 34 

2.6 Discussion of Results 35 

2.7 Conclusions and Recommendations 43 

HI. ANGULAR POSITION SENSOR ANALYSIS AND DESIGN 45 

3.1 Criteria Used in Selecting Operating Principle . . 45 

3.2 Position Sensor Description 45 

3.3 Discussion of the Dynamics of a Pendulum Augular 

Position Sensor 47 

3.4 The Influence of Acceleration on the Output 

Signal 51 

3. £> Transfer Funotion 54 

3.6 Switching Curve Computer ... 57 

IV. SIGNAL COMPARATOR 59 

4.1 General 59 

4.2 Experimental Analysis 60 

4.3 Discussion of Results and Conclusions. ...... 60 

V. DRIVING BISTABLE AMPLIFIERS WITH PROPORTIONAL 

AMPLIFIERS ' 63 

5.1 Introduction to Input Control Impedances of 

Bistable Amplifiers 63 

5.2 Problem Statement and Objective of Experiments * . 65 

5.3 Test Results end Conclusions 66 



iv 



VI. POWER AMPLIFIER AND BISTABLE REACTION JETS 

ANALYSIS ANL) DESIGN 69 

6.1 Determination of the Required Output Thrust .... 69 

6.2 Introductory Theory of Staging Unvented Bistable 

Amplifiers 70 

6.3 Initial Experimental Design 72 

6.4 Final Three Stage Amplifier Design 79 

6.5 Reaotion Jet Cascade Dynamios and Transfer Function 82 

VII. CONTROL SYSTEM SYNTHESIS, PERFORMANCE AND ANALYSIS. . . 86 

7.1 System Synthesis 86 

7.2 System Performance 91 

7.3 System Analysis 100 

APPENDIX 

A RATE SENSOR CALCULATIONS 103 

- A . 1 Rate Sensor Design Calculations 103 

A. 2 Calculation of First Stage Supply Pressure. . • . 104 

B CALCULATION OF THE REACTION JET TIMS DELAY CONSTANT . 106 

C SYSTEM ANALYSIS 110 

C.l Determining the Equation of the Switching Curve . 110 

C. 2 Limit Cycle Analysis • • . 112 

D PHOTOGRAPHS OF COMPONENTS 121 

BIBLIOGRAHTY 123 



v 



TABLE OF FIGURES 



Figure Page 

1.1 Phase Plane Representation of an Inertia Load 

Being Driven by a Constant Torque 4 

2.1 Vortex Angular Rate Sensor 6 

2.2 Sketch of the Vortex-Sink Angular Rate Sensor 9 

2.3 Cross Section of Probe 10 

2.4 Ideal Pressure Distribution Around a Circular 

Cylinder 12 

2.5 Sketch of Probe in Operation 14 

2.6 Effect of Viscosity on the Pressure Distribution 

Around a Circular Cylinder*. 19 

2.7 Typical Proportional Pure Fluid Amplifier 26 

2.8 Probe Boundary Pressure Versus No Load Flow 29 

2.9 Probe Output Pressure-Flow Characteristics 30 

2.10 Rate Sensor Steady State Output 36 

2.11 Rate Sensor Steady State Output 37 

2.12 Rate Sensor Frequency Response 38 

2.13 Typical Vortex Rate Sensor Frequency Response Data. . . 39 

3.1 Sketch of the Angular Position Sensor 46 

3.2 Analytical Model of the Position Sensor .48 

3.3 Sketch of the Position Sensor Mounted on the System . . 52 

3.4 Typical Steady State Position Sensor Output Showing 

Approximation to the Ideal Switching Curve 56 

4.1 Typical Summer-Amplif ier Test Results 61 

5.1 Bistable Fluid Jet Amplifier 64 

5.2 Input Characteristics of a Bistable Amplifier 64 

5.3 Effect of Biasing the Controls of a Bistable Amplifier 

on the Switching Pressures and Flows 67 

6.1 Control Volume for Determining the Output Thrust. ... 73 

6.2 Output Characteristics of the Two Stage Bistable 

Amplifier Without Driving Amplifier 77 

6.3 Output Characteristics of the Two Stage Bistable 

Amplifier ’With Driving Amplifier 78 

6.4 Reaction Jet Thrust Versus Supply Pressure 83 

6.5 Reaction Jet Input Characteristics 84 

7.1 Steady State Angular Rate Sensor Output While 

Loaded iirto the Summer-Amplifier • 88 

7.2 Block Diagram of Complete Control System 90 

7.3 System Cirouit Diagram 92 

7.4 Complete Control System and Demonstration Model .... 93 

7.5 System Response, Bang-Bang Control ;,.... 94 

7.6 System Response, Bang-Bang Control, Reduced Hysteresis. 96 

7.7 System Response, Bang-Bang Control, Without Rate 

Feedbaok 97 

7.8 System Response , Pulse Width Modulation Control .... 98 



vi 



Figure . Page 

7.9 System Response, P Vfl£ Control, Without Rate Feedback. « 99 

C.l Idealized Contactor Servomechanism Controlling 

an Inertia Load Ill 

C.2 Generalized Autonomous Control System Containing 

a Nonlinear Element Ill 

C.3 Polar Plot of G(jcJ) and -lAeq Without Rate 

Feedback 116 

C.4 Polar Plot of G(jcj) and -lA ea With Rate 

Feedback 116 

C.5 Bode Plots of Linear Transfer Funotion G(j60). .... 118 



D.l Vortex Rate Sensor Showing Details of 

Signal Amplification 121 

D.2 Summer-Amplifier 121 

D.3 Angular Position Sensor 122 

D.4 Power Amplifier Casoade and Reaotion Jets 122 



( 



vii 



CHAPTER I 



INTRODUCTION 



1.1 General Introduction to Pure Fluid .State Technology 

The basio concept of pure fluid amplification is that 
high energy fluid streams can be controlled by low energy 
streams. The concept was evolved at the Harry Diamond Lab- 
oratories of the U. S. Army Material Command, Washington, 

D. C. in 1959. * 

At the First Fluid Amplification Symposium held at 
the Laboratories in October of 1962, Brigadier General 
John G. Zierdt described the possible future of these compon 
ents : 



This advance has removed many of 
the limitations on conventional 
military and commercial pneumatic 
systems and has enabled an extreme- 
ly rapid extension of their capabil- 
ities. The use of low energy flow 
to direct and control high energy 
streams permits the design of pure 
pneumatic components, elements, and 
system^, capable of amplification, 
feedback, memory, logic functions, 
analogue computation, and digital 
computation. This, in turn, makes 
it possible to use;, pneumatic control 
systems, amplifiers, and computers 
where they could not be used before. 



1.2 rurpose ana Scope of Thesis , 

% 

Sinoe the inception of pure fluid amplification, much 



work has been done in analyzing individual components, re- 
cently the emphasis in the field has been turned to the de- 
velopment of complete systems. Although much work has been 
done along these lines, the quantitative data associated with 
such systems is scanty and inaccessible. 

One objective of this thesis is to develop techniques 
that may be. used in designing, analyzing and synthesizing these 
components, into a working control system. 

The specific objective is to develop a pneumatic con- 
trol system, using as few moving parts os possible, capable 
of maintaining a commanded angular position of a pure inertia 
load. Should errors be introduced, the control system is to 
reduce these errors to near zero in the shortest possible 
time within the limitations of the prime move used. 

Another design objective is to purposely make the 
system respond slowly so that the device may be used as an 
educational model demonstrating the working of individual 
components and the system as a whole. 

1.3 Discussion of the Proposed System to Attain the Design 
Objectives 

Minimum time control of a vehicle generally involves 
the use of th£ full force available from the prime mover in 
a bang-bang manner, i’or example, it is obvious that to move 
an automobile from a point, A, to another point, 3, in a min- 
imum amount of time, the operator should use maximum acceler- 
ation, at the right instant lock his brakes and slide into the 



2 



desired end point. The operator in this example would be 
using full driving force in one direction and full braking 
force in the other to achieve the minimum time objective. 

In controlling an angular inertia load, it con be 
shown that minimum time movements are obtained by using full 

1 - 4 * 

torque in one direction or the other . Figure 1.1 shows a 
phase plane portrait of an inertia load being driveh by a 
constant torque, whose sense is either plus or minus, for 
various initial conditions. 

It can be seen from Figure 1.1 that one curve exists, 
which, if used as a boundary to switch the sense of the torque, 

i ' 

will allow the system to reduce any errors to zero in mini- 
mum time. In order to switch the sense of the appliea torque 

i 

at the proper instant, it is obvi&us that the control system 
must be able to l) measure angular positions, 2} measure 
angular rate, 3) compute the coordinates of the switching 
curve, and 4) compare the values of angular rate and angular 
position with the values of the switching curve and switch 
the sense of the applied torque at the proper instant. 

The remaining chapters of this thesis deal with the 
design of the components needed to accomplish the above men- 
tioned functions and the synthesis of these components into 
a working system. 

* Superscripts Canute references in the Bibliography 
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Figure 1.1 

Phase Plane Representation of an Inertia Load, being 
Driven by a Constant < Torque 
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CHAPTER II 



VORTEX ANGUL..R RATE SENSOR ANALYSIS AND DESIGN 



2. 1 General Theory of Operation 

The principles on which a vortex angular rate sensor 
operate. are quite simple. Real difficulties quickly arise, 
however, when attempts are made at mathematical analysis of 
the complex flow which takes place in these devices or when 
attempts are made to obtain pressure or flow signals which are 
some function of the input angular rotation. It is primarily 
because of the latter difficulty that the methods of obtain- 
ing output signals have been, and still are, held confidential. 

The basic operation of a rate sensor is depicted in 
.Figures 2.]aend 2.1b. Figure 2.1a shows the flow state in the 

absence of angular rotation. In this. case a particle of 

% 

fluid passes from the supply plenum through the porous ring 
into the vortex chamber and out the exit hole. From mass 
continuity the radial flow velocity is described by 



U r 





2'rfr h 



where 



12 . 1 ) 



<4 = flow rate supplied to the sensor 
r = any arbitrary radius 
h s thickness of the vortex ohamber 
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Figure 2.1 

Vortex Aa pilar >ate Sensor 
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Thus just inside the porous ring, the radial velocity 



is given by 

u ° * as* 



and at the exit 



U* = 



JL 



2-rrR e h 



(2.3) 



where Q, and h are as given above 
and R e = radius of porous ring 
R 0 = exit hole radius 

Figure 2.1b shows the flow condition when the sensor 
is rotated at a constant angular rate. In this case a par- 
ticle of fluid entering the vortex chamber has the redial 
velocity given by equation (2.1) and a tangential component 
of velocity imparted to it while passing through the porous 
ring. In visualizing this tangential component, one must 
use a frame of reference other than the vortex chamber. Thus 
with respect to laboratory coordinates, the magnitude of the 
tangential component of velocity is given by 

U t0 = K 0 Gc> (2.4) 



where Uto = magnitude of the tangential velocity 
(■ d » magnitude of the input rotation 
Hence, a particle of fluid of mass dm will have angular mo- 
mentum given by Sq Uto dm. Conservation of Angular Momentum 
produces an expression for the tangential velocity lUt) at 
an- arbitrary radius (r): 

r Ut dm » R 0 Uto dm (2.5a) 

U t * R o 2 ^ (2.5b) 

r 
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For an ideal flow the familiar result of U t approaching infin- 
ity is seen as r approaches zero. 

Thus it can be seen that if a means can be devised to meas- 
ure the amount of tangential velocity directly or indirectly at 
a particular radius, this • velocity will change in proportion to 
the magnitude of CJ as to the direction of rotation. 

The method used to measure the smount of tangential veloc- 
ity present in the vortex chamber of the rate sensor used in this 
control system was developed by Sarpkaya 4 . The analytical analy- 
sis about to be discussed is also given by Sarpkaya but is included 

here to maintain continuity of this chapter and to set the back- 

♦ 

ground for the extension of Sarpkaya' s work which follows in the 
remaining articles*. , 

A schematic of the sensor developed by Sarpkaya is shown 
in Figure 2.2. This sensor differs from the one shown in Figure 
2.1 in that two sink-tubes are attached to the exit holes of the 
vortex chamber. These sink-tubes serve the following purposes: 

1) Tangential flow amplification - each tube can now be 
considered as a miniature elongated vortex chamber being fed from 
the top by a uniform radial flow on' which is superimposed the 
tangential swirQ. flow from the main vortex chamber, 

2) Directing radial flow - the. tubes receive the flow 
from the main vortex chamber and direct it toward the probe, 

3) The sink-tubes provide a housing for the probe. 

The probe, shown in Figure 2.3, consists of a small cir- 
cular cylinder divided into two chambers. Two small holes are 
drilled perpendicular to the center line of the probe and 



* To further credit the work of Sarpkaya, most of his 
symbols are maintained throughout this chapter. 
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Figure 2.2 

3k etch of -the Vortex- 3 ini: ‘inguiar Fete Sensor 
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Figure 2.3 

Cross Section of Probe 
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r 

aligned along a straight line drawn on the surface of the 
cylinder between its ends. The probe is inserted into the 
sink-tube perpendicular to the axis of flow, as shown in 
Figure 2.2. The holes are set at an angle to the oncoming 
flow of fluid in the sink-tube. 

The operation of this probe in sensing angular rota- 
tion is as follows: From potential flow theory the pressure 

distribution around a’ circular cylinder placed in a uniform 
flow field is given by 

P 0 + £ U p - . (1 - 4sin 2 ^J) ' (2.6) 

pressure in lbf/in 2 

static pressure in lbf/in 2 

density in lbf-sec 2 
in4 

axial sink velocity in in/sec 

angular position measured eround the 
cylinder from the stagnation point . 
n sketch of equation (2.6) is shown in Figure 2.4. In the 
absence of angulai* rotation of the device, the uniform radial 
velocity in the vortex chamber (U r ) turns into the sink-tube 
(now redesignated U s ) , strikes the probes and then passes out 
of the sensor. Because both holes are aligned at the same 
angular position along the probe, equal pressures exist in 
each chamber. The output signal of the probe, however, is 
the difference between these two pressures. Sinoe both 



?b = 
where = 

?o = 

P = 
9 -- 










Firrure 2.4 



Ideal Pressure J.istribution Around a Circular Cylinder 



12 



I 



chambers are at the same pressure, no output signal ( A 1) 
exists. 

No w let there exist angular rotation of the rate sensor. 
Figure 2.5 shows the velocity vectors in way of the probe. 

f 

The tangential velocity (Vp 0 ) now present in the sink-tube 
> 

has the bffeot at diverting the main axial sink flow. The 
top hole is thus placed in a uniform flow coming from an 
angle which is greater than that which existed without angu- 
lar rotation. The bottom hole is placed in a flow coming from 
an angle which is smaller. Thus, the pressure in the chamber, 
to which the top hole is connected, decreases, (point A in 
Figure 2.4), while the pressure in the bottom chamber increases, 
(point B in Figure 2.4). A pressure differential now exists 
across the probe which is a function of the magnitude of the 
angular rotation. It is readily visualized that had the angu- 
lar rotation been in the opposite direction that the reverse 
procedure would have taken place. Hence, the signal also 
senses direction of angular rotation as well as magnitude. 

It is assumed in this explanation that the relative magnitude 
of the sink-tube axial and tangential velocities are such 
that a long spiraling flow exists at the radius Tp 0 . 

To visualize the flow quantitatively and to obtain the 
function relating Ap ant COrefer again to Figure 2.5. By 
.inspection it can be seen that A Q, U 3 , and vpo are related 

. b y ; 



A$= tan" 1 - V PQ (2.7a) 

Us 
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Figure 2.5 

Sketch of Probe in Operation 

14 



if A0 is a small angle then: 

A0 - 122 (2.7b) 

u s 

* 

, To obtain the cazimun signal (A?) possible for a 
given CO , it is necessary to determine the maximum slope of 
the curve shown in Figure 2.4. Hence, differentiating 
equation (2.6) to obtain the slope gives: 



d P b _ P U 3 2 8 3in0Cos^ (2.8) 

3 0 2 

Now the derivative of equation (2.8) set to zero will give 
the 0 for maximum slope. 



g 

30 



2 




“ 8 P° s Cos 2 0 ,/9*'tt74 
2 

(2.9) 



Therefore, the probe holes should be located at a 
45 degree angle to the oncoming flow. This value substituted 
into equation (2.8) gives 




and in differential form for one hole: 
Ap| = 4 fpv a 




( 2 . 10 ) 



( 2 . 11 ) 



Equation (2.5b) gives the relation between the tangential 

l 

velocity existing in the sink-tube and the input angular 
rotation. With new symbols, equation (2.5b) beoomes: 



v 



po 



^Q 2 Cl) 
r po 



( 2 . 12 ) 
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Combining equations (2.7b), and (2.11), and (2.12) gives: 



A ? 

Pu s 2 



2 



4 R 0 2 CO 



U s^po 



for one hole (2.13a) 



Since two holes are being used equation ( 2. lX3a) becomes: 



A P 

P U 2 
' s 

2 
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po 



v> 

' (2.13b) 



This is the desired relationship between AP and CO in non- 
dimensional form. Note that for a given sensor and probe con- 
figuration, that E 0 and Tp 0 are fixed. Further, if the flow 
through the sensor is fixed, then U s is constant. Thus it is 
seen that equation (2.13b) is of the form: 



A ? = 


K O) 




($.14) 


where X - 


s p u 3 2 


h 0 2 w 


, a constant 
* depending on 
geometry and 
flow. 




2 


u s r po 



Data given by Sarpkaya bear out the linearity of this equation 
over a substantial range ( of CO . As might well be expected, 
from inspection of Figure 2.4, linearity is preserved provided 
A Q is small. 

Thus far the equations derived to describe the opera- 
tion of this rate sensor, are based on an idealized potential 
flow. It is now appropriate to consider the effects of vis- 
cosity on the operation of the sensor and the describing equa- 
tions. Two major effects of viscosity will be noticed: 
l) Viscosity will cause the vortex motion in the 
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vortex chamber and sink-tube to deviate from the ideal vortex 



given jy equation (2.5b). 

2) Viscosity will cause the flow around the probe to 
separate resulting in an alteration of the pressure distri- 

v 

but ion around the proqe given by equation (2.6). 

Effect 1) is studied in great detail by Sarpkaya, by analyzing 
the boundary layer equations in the chamber for the radial flow 
with small tangential flows superimposed on the main flow. 

The results show that viscosity reduces the strength of circula- 
tion as the flow travels from R 0 to rpo« Sarpkaya introduces 
a factor he colls the viscous efficiency of the sensor defined 
as: 




where = strength of vortex circulation at 
radius r 

P 0 = strength of vortex circulation at 
radius R Q 

The viscous efficiency between the porous ring and the probe 
is: 

r 

E — — ~~ - v no r po - v po r PQ 

* o R Q U 2tT R q 

(2.16a) 

from which 

V - _ R b ^ E (2.16b) 

r po 
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Wj.th equation (2.16b) used in lieu of equation (2.12) and then 
combined with eouations (2.7b) and (2.11), the new relation- 
ship between AP and CO becomes 




8 for two holes. (2. 

U a r 



Intuitively one might expect thi.s efficiency to be relatively 
high provided the circulation is quite small. This is shown to 
be true by Sarpkaya in a plot of E versus Reynolds Number. 

Effect 2) is an extremely important consideration which 
seems to have been neglected by Sarpkaya in his analysis of 
the sensor. Non ideal flow around the probe radically re- 
duces the slope of A P versus Q in the vicinity of 45°. 

For the case of ideal flow around the probe, the nondi- 
mensional slope is given by equation (2.11). 

2 >' ■ 



ap|/ Bh 



= ^ 






( 2 . 11 ) 

repeated 



From actual measurements taken on the probe used in 

f 

# I 

the sensor constructed by the author, this slopb was found to 
vary between 2.58 and 3.17 for the flows considered. The ef- 
fect of viscosity is more clearly shown in Figure 2.6 where 
nondimensional pressure distribution for two actual flows 
versus 0 is shown along with that predicted by potential 
flow theory. Thus the relationship between AP and CO with 
the introduction of the new slope of Ap vs A $ becomes : 
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/°Us/2 | 




Firure 2.6 

Effect of Viscosity on the Pressure Distribution 
Around a Circular Cylinder 



19 



Ia. 



/0U 3 2/2 



5.16 
u s r po 



for two holes. 

(2.18) 



2.2 Rate Sensor Design and Design Procedure * 

The author did not consider the effect of slope reduc- 
tion due to viscosity of the Ap versus Q curve, discussed 
at the end of article 2.1, until after the rate sensor had 
been designed. Hence, the equations used for the design were 



Ap = 
A0 = 



8 



* / 
R 0 2 0)2 



r po U s 



A$> two .holes (2.19) 

( 2 . 20 ) 



Prom the design specifications discussed in the intro- 
ductory chapter to this thesis, O^qx was selected to be 
70 degrees per second. In order that the output of the sen- 
sor be linear over the range - Q) rflT 6Q) £0-L ax > tl10 value 
A Q must be constrained to stay within the linear portion 
of the AP versus Q curve. Analysis of the data given by 
Sarpkaya has shown that A Q should not exceed roughly ±3° 
for linearity. <' 

Inspection of equation (2.20) shows that R 0 , r 0Q> E, 
and U s are still to be determined. R Q was chosen on size 
considerations of the sensor and was selected to be three 
inches. The value of U s will depend on the magnitude of the 






See Ajppenaix A.l for detailed numerical calculations 
in this section. 
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desired output signal AP and can be determined from equa- 
tion (2.19). Once U s is determined, ^ can be calculated from 
the relationship 



4/2 - 



'TT'r Q 2 - 2c(2r q ) 



Uc 



( 2 . 21 ) 



where r< 



radius of sink 



c = radius of probe 



In order to limit 4 to a reasonable value, it is obvious that 
the magnitude of AP must be kept relatively small. For this 
design a value of Ap = 0.1 psi to occur at wqs 

selected. 

v 

Having determined 4» the Reynolds Number can be deter- 
mined from 



He 



where y 
h 



A ft/ 3 ... (2.22) 

V V 



J± 

p 



kinematic viscosity 



thickness of vortex chamber 



With the value of Re determined, the experimental data given 
by Sarpkaya on viscous efficiencies can be entered and a value 
of S determined. 

Knowing all the values of the variables in equation 
(2.20) except r 00 , allows one tor solve for this variable. The 
position of the probe holes from the center of the sink-tube 
(distance a in Figure 2.3) can then be determined from the 
relationship 



21 



( 2 . 8 3 ) 




which can be verified by inspection of Figures 2.3 end 2.5. 
One word of caution must be made concerning a minimum value 
of r n0 . Viscosity plays an increasingly important role as 
r approaches zero in a real vortex. If r 00 is too large, 
full advantage is not made of the amplifying properties of 
the vortex. On the other hand, if it is made too 3mall, vis- 
cosity rapidly decreases the strength of the vortex. This 
phenomena' remains to be investigated more fully and is men- 
tioned here only for design considerations. 

The design data used in constructing the rate sensor 
is summarized in Table 1. The material used in construction 
was plexiglass. The porous ring was mace of felt faced with 
fine wire screening which gave rigidity to the ring and acted 
as a flow straightener . The probe was constructed by cutting 
a 1 3/4 inch piece of number 19 gauge stainless steel hyp'o- 
dermic tubing in half and inserting a plug. The number 80 
Drill Size holes were then drilled 0.095 inches apart. This 
construction allowed the holes to be aligned by rotating one 
side of the probe relative to the other. 

2.3 Signal amplification 

t 

The magnitude of the signal output (pressure and flow) 
by itself is too small to be used conveniently in the control 
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TABLE I 



w 

SUMMARY OP RATE SENSOR DESIGN DATA 



Specified Data: 



k-) max 


70°/sec 


Ho = 


3 in. 


^ P max “ ’ 


0.1 lbf/in 2 



a 9 ^-' 30 

Arbitrarily Chosen Data: 

•r 3 = 0.1562 in. 

c - 0.0212 in. (//19 gauge hypodermic 

tubing OD x .ID = 

0.0425 x 0.0270 in) 

probe hole diameter = 0.0135 in. (rfQO Drill) 

f 

Derived Data: 

% 

r * 0.0512 in. 

a r 0.0475 in. 

U 3 = 171.2 ft/seo 

Q - 9.05 ft 3 /min 

S = 0.50 



I 
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system loop and therefore must be amplified. The state of 
the art of pure fluid signal amplification, though not yet 
mature, is advanced to the point where a combination of ex- 
perimental and analytical methods can be used to obtain rea- 
sonable signal amplification with tolerable signal to noise 
ratios. 

The design problems faced at this stage of the develop- . 
-ment of the rate sensor ere: 

1) Given the output pressure-flow relations of the 
rate sensor probe, whet are the supply pressure settings 
needed at each amplifier stage to insure linearity of the 
signal? 

2) How many stages of amplification are needed to 
raise the pressure and f low- to acceptable levels? 

The last question could not be answered at this 
stage of the design for the rest of the systems com- 
ponents had not yet been designed. Therefore, it was 
decided to use two stages of amplification and after 
the amplifier stages were attached, to accept the out- 
put pressure-flow relations of the final amplification 
stage as the output characteristics of the sensor. 

The first question, although a relatively straight for- 
ward design question, is covered only lightly in the literature. 



The procedure used to obtain the supply pressure settings 
described below is by no means exact but it will let the 
designer set the pressures close to their optimum settings 
on the first try. 

The amplifiers used to amplify the rate signal output 
were Corning Class Work* s -Standard Size Proportional Amplifiers. 
The important dimensions of these devices are: 

power nozzle width a w s 0.010 in. 
control nozzle width = 1.5w 

depth - 2.5w 

» receiver distance =• lOw 

Figure 2.7 is a sketch of the interaction region of a typi- 
cal proportional amplifier which operates prinoiply by mo- 
mentum interaction. Following the development given by 

5 

bexter , the deflection of the power stream due to one control 
momentum only is determined as follows: 



\ 

t 




where r m c V Q 

m c +mp 

V a = 

mc + m p 

and m p « mass flow rate of the power jet 

I 

i mass flow rate of the control jet 
V c s velocity of the control" jet 
Vp i velocity of the power jet 
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. *. tan Q ~ . 

but mv = PaV 
tan 



m c 



a 



■P Vr 



- P*c V 

pA n V ‘ 

' P P 



Pdv; c V 0 2 
P d w p V p 2 



where d r 
w c = 
w p = 

tan$ - 



depth 

control nozzle width 
power nozzle width 



w c V 



w_ v 2 
P v p 



but V 2 = 2/p 

/. tan $ - . 



w. 



2 /P Pc 



» P 2/ P ? i 



V?, 



w p p p 



(2.24) 



(2*. 25) 



Reilly and lioynihan p ® indicate that the output of a 
proportional device w ill remain linear provided the deflection 
of the power jet at the receiver ports does not exceed one 
half the nozzle width. Thus: 



tan Q ^ **?/% 



for linearity 



where L = receiver distance 
Combining this with equation (2.25) gives: 



w. 



2L 



w c P C 

w p p p 



for linearity 
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This result is readily extended to two control ports with the 
result: 



v ‘p y ;v c ( r cL-^"cR ) for 
2L ” v/p linearity 



Having developed the basic equations of the amplifier, 
the output characteristics of the rete sensor probe nmst now 
be investigated. According to potential flov; theory, if the 
probe is set so that the probe holes are at 45 degree angles to 
the oncoming flow, the pressure in the probe chambers is. below 
ambient. Since the probe chambers are connected to the control 
ports of the first stage amplifier, this means that only the 
pressure differential existing in the two chambers and com- 
municated to the first stbge control ports is available for 
deflecting the jet. This is undesirable as it results in 
less gain when compared to operation with positive flow through 
both control ports. However, due to viscosity the ideal flow 
pattern is radically altered so that the pressure existing in 
the chambers when the flow angle is set to 45 degrees is posi- 
tive. 

The various flow rates to be expected out of either 
end of the probe, as the main sink flow is deflected, can be 
approximated by setting the sink flow at its operating value 
and turning the probe with respect to the flow. The results 
of such an experiment are shown in figure 2.8. 

The output pressure-f low characteristics for a given 
are shown in Figure 2.9. An infinite number of such curves 
exist since varies with A Q* 
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Figure 2.3 

Probe boundary Pressure versus I’o Load Flow 
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4.0 

of water) 



-figure 2.9 

Probe Outnut Pressure*- Plow Characteristics 



The procedure initially used to set the first stage 
supply pressure (?^) was to determine the output pressure- 
flow relationship for maximum A ? and superimpose the input 
characteristics of the first stage amplifier to determine 



(1 



cL 



"cR^max 



nax . Then, from equation (2.25), the 
supply pressure was determined. However, when this procedure 
was followed, it was found that, because the control port area 
of the first stage amplifier was so large compared to the probe 
hole area (nearly three times as large) , the input impedance 
of the amplifier was nearly zero. 

Circumventing this problem, equation (2.24) was used 
in lieu of (2.25). With the constraint placed on tan Q and 
for two control port operation, equation (2.24) becomes: 



v; 



p 



, w 0 (V cL 2 - V cS 2 ) 



max 



(2.27) 



2L 



w p V d 2 



From measurements of the flow rate occurring at the 

expected maximum pressure just outside the probe hole, 

(V ^ 2 - V ^ 2 ) was obtained from V rtT 2 - VAj* 2 c 1 
v cL cR max Ci - 0iX £ 1 



(Q C L 2 aex - Q. c tj 2 m i n ) which then allowed Vp to be obtained from 

i 

equation (2.27) 

To obtain Pp the following relationship was used: 



<i T 



•where Q. 



'dp 



A P V P 



\fz7p p. 



(2.28) 



'dp 



discharge coeficient of the power 
nozzle which was determined from the 
design data sheet published by Corning 
glass 'Works. 
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From the calculations in Appendix A.?,, was determined to 
be 7.7 inches of water. 

The second stage pressure was set using information 

r 

given in the design data sheet published by Corning Glass 
'.York's for these amplifiers. This sheet states that for max- 
imum gain the control bias level pressure should be one tenth 
the power jet pressure*. 

The bias level of the second stage controls was deter- 
mined by setting the first stage supply pressure at 7.7 
inches of water and measuring the pressure in the line connect- 
ing the output of the first stage and the control port of the 

m 

first stage. This procedure resulted in an initial setting of 
Pp 2 equal to 5 psi . 

2.4 Hate Sensor Dynamics 

This rate sensor is characterized dynamically by 1) 
a pure time delay resulting from the time it takes for the flu- 
id particles containing tangential momentum to reach the probe, 
2) a first order lag due to the resistances of the probe holes 
and control nozzles of the first stage amplifier and the vol- 
ume connecting them, and 3) the dynamics of the staging am- • 
plifiers and the lines used in interconnection. 

The transport delay is the most significant. It may 

$ 

* This procedure used in setting the first stage 
pressure resulted in a Pp-^ » 1.95 inches of Vrnter and hence 

was not U3ed. 
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be approximated by the time it takes for the sensor to com- 
pletely renew the supply of fluid in the vortex chamber^ 
and sink-tubes. Thus: 

7* - — — seconds (2.29) 

where V - total volume 

Q, = flow rate given by equation (2.23) 

The total volume (V) is given by: 

V - 'TfH 0 2 h •+* 2'n'r s lp 

where x? Q , h and r s are as before 

1 r distance to probe measured along the 
^ axis of the sink-tube 

For the geometry and flow rates of this device, the value of 
T is 0.0556 seconds. 

The time constant of the first order lag referred to 
in item 2) was found to be extremely small and was neglected 
in arriving at the transfer function of the rate sensor. 

The dynamics of the staging amplifiers are discussed by 
Belsterling and Tsui^ and Boothe?. These works show the de- 
vices to respond very much faster than the time delay of this 
rate sensor. In addition, the design data sheet for these 
amplifiers indicate that the frequency gain is over 1000 cps, 
well out of the range’ of the rate sensor. 

Thus the transfer function relating the output (A?) 
to the input angular rate [CO) is given by: 
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t 



A? 




Oj 



where • 



steady state gain 



7 % 



pure tine delay given by 
equation (29) 



Laplace operator d/dt 



2 . 5 Testing Procedure and Results 

The rate sensor was made into a self _ conta ined pack- 
age by attaching the staging amplifiers directly to one side 
of the vortex chamber. The sensor was equipped with a man- 
ifold allowing air to be supplied with only one line. This 
manifold pressure' was then reduced to the required pressures 
needed at the vortex chamber supply plenum and power jets 
of the staging amplifiers. Variable resistors were used to 
set the pressures and were constructed from tapered dowel 
pins through which holes were drilled. Photographs of the 



assembled rate sensor are shown in Figures 7.4 ana D.l . A 
differential pressure to electric transducer was attached to 
the outputs of the second stage amplifier and the assembled 
device was tested on a rate table at the- Instrumentation 
Laboratory of the Massachusetts institute of Technology. 



1') Steady state test over the range 1.10— 60 ^1.15 
radians per second with the output loaded by orifices approx- 
imating the size of the control parts of an amplifier similiar 
to those used on the rate sensor for signal amplification. 



The testing procedure consisted of the following: 
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2) Steady stete-test over the sene range with the out- 

* 

put flow equal to zero; i.e., infinite load. 

3) frequency response test at discrete frequencies 
up to 4.67 cps. 

4) Steady state test while the device was being simut- 

I 

taneously rotated about its sensing axis and an arbitrary axis 
The test results for the first three tests are shown in 
Figures 2.10, 2.11, 2.12, and 2.13. 

2.6 Discussion of Results 

The results of the steady state rate tests show the 
following: 

1. The very snail signals available fron the probe can 
indeed be amplified into reasonable signals. 

2. The output of the sensor is not quite linear over — 
; the range tested. 

w 3. The gain of the sensor with clockwise rotation is 

* significantly larger than the gain with counter clockwise ro-* 
tat ion when the bias is subtracted out of the results. 

4. The gain increases continually as the input rate 
is increased. 

5. The dog leg effect produces a highly undesirable 
rapid change of gain in the counter clockwise direction. 

In an atte&pt to' explain the last four items, accurate 
flow and pressure measurements were taken on each side of 
the probe as the angle of incidence of the probe holes with 
the oncoming flow was varied. Originally, in aligning the 
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Rate Sensor Steady State Output 
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Sweep = 0.2 sec /cm 
Amplitude Ratio = -1.4l db 
Phase Shift = 62° 
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Frequency >= 4.67 cps 
Sweep = 0.1 sec/cm 
amplitude r Ratio = -3.93 db 
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Input Signal (<^ nax = 1.02 rad/sec) 



Figure 2.13 

Typical Vortex Rate Sensor Frequency Response Data 



39 



I , 

probe holes (the turning of ono probe hole with the other held 
fixed), only the dead-ended pressures were used, taking the 
probe holes as being aligned when the pressures were equal. 

As v/as pointed out in article 2.3, signal amplifica- 
tion is essentially dependent on the flows coming from the 
probes. When these flows were measured, it was found that 
equal dead-ended probe pressures resulted in radically dif- 
ferent flows indicating that, although the holes were sup- 
posedly of the same area, they were behaving as if they were 
not. 

The probe was removed from the sensor and observed under 
a 50 power microsoope. This showed that during the machining 
process, the drill had been allowed to cut a trough into the 
entrance to one of the holes. The effect of this trough was 
apparently to catch more flow from the main stream resulting 
in different orifice characteristics. 

The effect of having the two probe holes producing 
different flows for equal boundary pressures was as follows: 

It has previously been shown that, as the sink-tube axial 
velocity is deflected, the pressure outside the probe holes 
increases or decreases depending on the sense of angular ro- 
tation and the probe hole in question. As the pressure out- 
side the probe hole which gives less flow for a given pressure 
is reduced, a point is reached where the flow approaches zero. 
By referring to Figure 2.8, it can be seen that, as the flow 
approaches zero, the relation between the pressure outside 
the hole and the resulting flow beoomes highly non-linear. 
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This highly non-linear flow region produces the dogleg ef- 
fect noted above. This effect is not present when the direc- 
tion of rotation is reversed for, as the boundary pressure 
outside the hole which gives the greater flow for a given 
pressure is reduced, the non-linear portion of the curve is ■ 
avoided*. At the same time, the other probe hole is experi- 
encing increasing pressure and hence this combination produces 
no dogleg effect. Thus, to remove the effect, the probe hole 
without the trough was altered to give approximately the same 
pressure-flow characteristics as the other probe hole. 

To explain the difference in the maximum signal output 

between clockwise and cqunter clockwise rotation, one has only 

to realize that the positioning of the probe holes with respect 

to the oenter 'line of the sink-tube is critical. If the holes 
* 

are off center the slightest amount, the output will be de- 
creased on one side and increased on the other resulting in 
different peak outputs. 

The effect of the gain continually increasing with in- 
creasing CO can be partly explained by realizing that deflec- 
tion of the first stage power jet is roughly proportional to 
( A~i) 2 . 'Generally the output pressure of a proportional am- 
plifier is given by 

Ap 0 s kjAp 0 

but Ap 0 = ic 2 (A >i 2 ) 

* At least over the range of rotations tested. 









hence A ? 0 * (A«c ) 2 

^ » 

Now as shown in Figure 2.8, A *4 is approximately proportional 

t 

to CO over a substantial range of operation, making the output 
pressure generally proportional to the square of CO • This 
of course is only approximate but does roughly explain the 
increasing gain of. the device. In fact, if the output pres- 

f 

sure is converged to flow using the output characteristics 
of the second stage amplifier , one finds the output flow quite 
linear with respect to co . 

Although accurate verifications of the explanations 

1 

given for items 2, 3, 4, and 5 above were not accomplished 

due to inaccessibility of the rate table, the explanations 

g^ven were partially verified and showed that the dogleg had 

been removed but that only slight improvement was made in 

balancing the outputs by centering the probe. 

The zero input bias signal was found to be only partly 

due to the unequal flows from the probe. The amplifiers 

* 

themselves are biased when the control flows are zero. This 
can be reduced somewhat by proper arrangement of the cascade 
so that the output bias of one tends to cancel the output 
bias of the next stage. However, adjusting orifices general- 
ly must be placed somewhere in the cascade to remove the bias. 

The results of the frequency response tests verified the 
pure time delay nature of the dynamics of this sensor. The 
experimental phase data is approximate at the low frequencies 
because the method used to obtain the data (the input and 

A 

V2 



output were recorded on permanent recording paper and were 
compared in order- to measure phase) could not discern small 
phase shifts. The amplitude data, however, is quite accurate. 

Lissajous figures on an osciloscope could not be used 
to obtain frequency response* data because the output noise 
moved the figures around the scope precluding any measurements*. 

The best data was obtained by taking pictures of the 
oscilloscope traces and measuring phase ana attentuation data 
from these. The rapid increase in phase verifies the existence 
of a pure time delay. The deviation from theory in phase 
shift and attenuation is probably due to volumes in the lines 
connecting the first and second stage amplifiers and the vol- 
ume and resistance of the orifices between the output of the 
second stage amplifier and the pressure transducer. 

The results of the last test showed the output signal 

f 

did not change while the sensor was being simultaneously 

i , 

rotated about its sensing axis and an arbitrary axis. 

2.7 Conclusions and Recommendations 

Throughout the development of this rate sensor little 
attempt has been made at optimization. The objective was to 
develop a rate sensor capable of being used in a control system. 
To this end the following conclusions are drawn: 



* The nol'se of fluidic proportional amplifiers is of 
a peculiar nature. Although high frequency noise is prevalent 
as shown in Figure 2.13, a more disturbing very low frequency 
noise is present which causes the signal to wander and make 
unexpected excursions which appear almost as steady state bias 
changes. 
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1. Use of the design procedure suggested in article 2.2 
will allow the designer to theoretically design a rate sensor 
peculiar to his applications. 

2. The steady state gain of the unamplified signal 
may be predicted analytically but due to the low input im- 
pedance of the type of amplifier used, the amplified gain will 
have to be determined experimentally. 

3. The dynamics of the sensor may be determined to a 
reasonable degree of accuracy by theory. 

4. The output of the sensor is very dependent on the 
accuracy of the machining of the probe holes and alignment in 
the sink-tube. 

5. Noise of the very low frequency type is a serious 
problem causing unexpected wandering of the DC level of the 
signals . 

6. Contrary to the results of Sarpkaya, it was found 
that closing one of the sink-tubes and keeping the flow con- 
stant, resulted in an increased output signal for a given CO . 
Since this equivalently means nearly the same signal with less 
power consumption, the author recommends that sensors of this 
type be built with only one sink-tube and, as recommended by 
Sarpkaya, that the end opposite the sink-tube be designed to 
streamline the flow into the sink-tube. 



CHAPTER III 



ANGULAR POSITION SENSOR ANALYSIS AND DESIGN 



3.1 Criteria Used in Selecting Operating Principle 



Angular position of the system can be obtained by 

i 

analog 8 or digital 3 ’ 9 integration of the angular rate signal 
available from the vortex rate sensor. Digital integration 
would possibly involve developing a pressure sensitive os- 
cillator and some form of digital-to-enalog converter, both 
of which are beyond the time limitations of this thesis. 

Since a linear position signal over a large range of angular 
positions for demonstration purposes was desired, analog 

r . 

integration was ruled out since very large 'volumes would have 

« 

been required to"s tore 1 ' the position sighal. furthermore, the 

7 

state of the art of pure fluid analog computation is such 

that up-down integration is inaccurate. ; 

♦ 

. t 

3.2 Position Sensor Description 



i 

Because of the problems mentioned above, it was decided 

■i 

to use a moving part position sensor consisting of a pendulum 
set between a pair of flapper nozzle valves to obtain a dif- 
ferential pressure signal that is a function of angular posi- 
tion. A schematic of this sensor is shown in Figure 3.1 and 
a photograph of the final design is shown in Figure D.3. The 
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x’icure 3.1 

o’retck c" the ^nrulor i os it ion Sensor 
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device is designed so that as the sensor rotates, the pendulum 
trys to return to a vertical position. In so doing, it widens 
the gap on one side anc closes the gap on the other until the 
moments caused by the flow forces from the nozzles, balance 
the moment due to the weight of the pendulum. .Vhen this 
occurs a differential pressure exists in the chambers and 
constitutes the output signal of the sensor. 

It is readily seen that the output differential pres- 
sure signal varies as the sine of the angular position and 
that the sensor also senses direction of angular position. 

3.3 Discussion of the Dynamics of a Pendulum Angular 
Position Sensor 

If the bearing which supports the pendulum is relative- 
ly free of damping, the angular position of the pendulum is 
unstable and the slightest perturbation will start the device 
oscillating. In the bang-bang control envisioned for the 
final system, this oscillation is undesirable. To understand 
the reasons for this instability, a linear model of the system 
can be constructed to describe the behavior of the system 

f 

when the angular position of the pendulum is small. Figure 3.2 
shows an analytical model of the sensor. Summing moments 
about the pivot point gives the following differential equa- 
tion: 




1 te* bl b $* lia n + Vl) & 



T d U) 



4 



\ 



47 



///////// 



! 




v iRure 3.2 

Analytical I'odel of the Focition densor 



J+8 



where 



9 ‘ 

9 = S 2 ^/dt 2 

bl^ = equivalent viscous damping 

k - equivalent spring constant of the 

flapper nozzle valves 

T^ = a disturbance torque 



Neglecting the term \lQ , since it is much smaller than 



and introducing the symbols J* 
X= kl n , and s=~d/dt gives • 

Js 2 + E.s + X Q = 

• m 

The spring force X ^is given by 



V/ 2 

T w ’ 



B= bl 



b» 



(3.1J 



s (P X A - P 2 A)l n = A?Al a 

where P]_ = pressure in one chamber 



pressure in the opposite chamber 



rearranging gives 



X 



Al A? 
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(3.2) 



The linearized transfer function between A? and$ 
is, in general, a first order leg of the form‘d. 



Ap 

B 



(3.3) 



Ts+1 



^9 



V 



Substituting equation (3.3) into (3.2) gives 

k r A1 n°l = C 2 . (3.4) 

T S4-1 Ts+1 

* 

If this equation is substituted into equation (3.1) the 



result is 



[ 



Js 2 + 3s +2L 



e = 



Tsfl J 

From this the transfer function between Q and is 



$ = Ts ± 1 

T d JTs 3 +(J+BT)s2 f BSfC 2 

The characteristic equation of this transfer function is 

j^JTs 3 + (J + BT) s 2 + £3 + C 2 ] 

The South stability criterion for linear systems gives the 
condition for stability as 



B(J*BT) > JTC 2 



or 

* 

«• 

B J 2 + 4JT " J (3.5) 

— 

where the positive sign for the radical has been chosen since 

0 

the damping can not be negative. If equation (3.5) is re- 
written as 



B -> 1/2 






f 4JC 2 






r 
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it can be seen that the value of Cg is the important para- 
meter determining the stability of the system. As C 2 in- 
creases, 3 must be made correspondingly larger. From 
equation (3.4), Cg is the value of the steady state spring 
constant of the flapper-nozzle valves. Thus, as this spring 
constant increases, the magnitude of the damping must be 
increased. 

Initially, to introduce damping into the system, the 

ball bearing used to support the pendulum was packed with 

heavy grease. This did not introduce a sufficient amount 

* « 

of damping into the system. The problem was finally solved 

by placing two small squeeze film dampers on either side of 

the pendulum and adding a minute amount of grease to aid in 

* 

the damping process . 

3.4 The Influence of Acceleration on the Output Signal 

The output signal of this sensor will be influenced 
by. the angular acceleration of the device to which it is 
attached. Figure 3.3 is a sketch showing the position sen- 
sor mounted on the system and offset from the axis of rotation 
by a distance r. 

If the sensor is given an angular acceleration & , 



* Over most pressure settings the air in the squeeze 
film dampers provided sufficient damping. However, if the 
chamber pressures were made too high then the size dampers 
used were inadequate. 
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then the end of the pendulum receives a linear acceleration 
rCL. The torque on the pendulum is given by 

T - JH ra (3.6) 

g 

* 

This torque must be balanced by the pressure forces on the 
end of the pendulum giving 



Ap Jji A1 n : -g- r a. 

t 

or 

Ap la = -as a l3 - 7) 

This A? corresponds to an error signal. It may be reduced 
by making r small assuming the other parameters are fixed. 

The moment of inertia (I) of the pendulum is given by 



I 




Thus the resulting torque on the pendulum is given by 



T AA z -t V s a (3.8) 

This torque too must be balanced by the pressure forces on the 
end of the pendulum giving 

^ ^AA g CL 



or 




wl 8 w 

gAl n 



a 



(3.9 ) ■ 



* "Subscripts LA and Aa refer 'to linear acceleration 
and angular acceleration respectively. 
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This A P also represents an errot in measuring angular position. 
The total error signal is the sum of the two differential pres- 
sures given by equations (3.7) and 3.9) 



Ap 

.error V' 



J wr 

r] ^i n 



f wl 2 W 
gAl n 



a (3.10) 



This pressure error may be converted into an equivalent angu- 
lar position error by the equation 



d z Sin" 1 I A p (3.11) 

\ wl w / 

l 

which is obtained by a steady state torque balance on the 
pendulum when at an angle^ from the vertical. Substitution 
of the value of A? from equation (3.10) into equation (3.11) 
will give the error in terms of angular position. Primarily 
because the accelerations to be cea&t with in this system are 
quite small, substitution of the actual values into equations 
(3.10) ana (3.11) showed that the maximum probable error would 
be less than one degree. If this system were to be made into 
a very fast high torque system, this type of position sensor 
would probably be undesirable unless the designer wished to 
have acceleration feed back in the system. 



3.5 Transfer Function 

Because the amount of damping added to the system to 
xaake it stable is difficult, to determine, it was decided to 
connect the output signal into a pressure to electric trans- 
ducer anc experimentally determine the transfer function. By 

r 

observing step responses on an oscilloscope, the system dynamics 



can readily be determined. 

The procedure was as follows: 

v 

1) The output signal was connected to the transducer 
and loaded into orifices which approximated the size to be 
expected at the control ports of the summing device. 

2) The output signal was then displayed on an oscil- 
loscope. 

3) The pendulum was given an initial deflection and 
then released. 

The response of the pendulum was observed to be es- 
sentially that of a first order lsg. 3y measuring the time it 

I 

took the signal to decay to 36.8 percent of Its initial value, 

I 

gave a close approximation to the time constant of the sensor. 
3y averaging a large number of such measurements the time con- 
stant of the system was found to be 0.04 seconds. This time 
constant was found to be much larger than the time constant 
associated with the flapper-nozzle chambers*. Apparently, 
the measured time constant is essentially that of the mass- 
damping time constant of the system. 

A typical steady state output curve is shown in 
figure 3.4 from which the steady state gain may be measured. 

The’ transfer function, relating the output pressure 
signal to the input angular position signal for small valves, 
of Q , is given by 



* The time constant cue to the volume of the flapper 
nozzle chambers and connecting lines was found to be 0.008 
seconds. The calculations are not shown here as they can be 
found- in any standard reference such as reference 10. 
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'■i'-ure 3. + 

Epical Steady State Position Censor Outnvt 
Showinc Approximation to the I^cal itchinr, Curve 
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A? 

6 



(3.12) 



r K 2 
r 2 a *l 

where K 2 - gain to be determined from a plot 

similar to Figure 3.4 once the final 
system pressures have been set 

7*2 = 0.04 seconds 

Before concluding this chapter it is to be noticed 
that little information has been given concerning design 
procedures. This was purposely done as such procedures are 
more adequately covered in the literature 1 ^. However, the 
designer must keep in mind the over all system and the re- 
lationship of each component to the others. For example, 
it is obvious that the output of the position signal is to 
eventually be compared with the rate signal. It is therefore 
reasonable to assume that these signals should be of the 
same order of magnitude. 

3.6 Switching Curve Computer 

I 

As shown in Appendix C.l, the ideal switching curve 
for this system is a parabola passing through the origin. 

The slope of this curve is infinite at the origin. Since it 
was planned to generate the switching curve by operating on 
the output signal of the position sensor, this meant that as 
0 varied slightly from its zero position an infinite gain 
would be required. Because of this consideration and the 
fact that the output signal of the rate sensor is not depend- 
able at very low angular rates, the output of the position 

h 
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sensor, without modification, was used as the switching curve. 

Such a sine curve output can be adjusted to approximate 
a parabola over a large range of position signals as shown in 
Figure 3.4. The approximation is poor when Q is small, how- 
ever, and thus, it is to be expected that "linear switching" 
will take plaoe when the errors introduced into the system 
are snail. 
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CHAPTER IV 



SIGNAL COMPARATOR 



4.1 General 

The method^ selected to algebraically sum the pressure 
signals from the rate sensor and position sensor consists of 
adding additional sets of controls to a proportional amplifier 
and having the power jet deflect in accordance with the vector 
sum of their momenturas*. A stream interaction proportional 
■amplifier, operating on this principle, is available from 
Corning Glass Works** and has two sets of control ports. The 
design data sheet accompanying this amplifier does not in- 
dicate that 'the second set of controls can be used for sum- 
ming and, therefore, gives no data on the gain of these ports. 
In fact, the ports are not equipped with fittings and, there- 
fore, had to be adapted for use in this oontrol system. A 
photograph of this "summer-amplifier" is shown in Figure i>.2. 
The vertical controls are the secondary controls and were-' 
adapted for use by inserting hypodermic tubing into the holes 
and then cementing the tubing into place. 



* Other methods of summing signals, using no moving 
parts, are available and could have been used in this system. 

** Catalog Number FE2511-2-1211 September, 1965 
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4.2 Experimental Analysis 

A silhouette of the amplifier described above is given 
in the design data sheet. From, this silhouette, it v/as found 
that the width of the secondary control ports is less than 

4 

the width of the primary ports. In eddition, the secondary 
ports are not perpendicular to the power jet. Because of these 
geometrical differences, it is obvious that if the same pres- 
sure differential is applied to both sets of controls ports, 
the power jet will be deflected a greater amount by the pri- 
mary control ports. As a consequence of this greater deflec- 
tion, the output pressure will be larger. Therefore, the 
primary control ports will exhibit a greater gain. 

To determine, the difference in the gain of these two 
sets of control nozzles, a number of tests were performed on 
the amplifier. The tests consisted of 1) systematic varia- 
tion of the control port bias levels about the levels ex- 
pected from the outputs of the rate and position sensors, 

2) systematic variation of the supply pressure, and 3) at any 
particular bias level and supply pressure setting, recording 
the gains of the two sets of control ports. 

Typical data from these tests' for one particular com- 
bination of supply pressure and bias level settings is shown 
in Figure 4.1 During all these tests, the output of the 
summer-amplifier was loaded into a manometer. 

4.3 Discussion of Results and Conclusions 

The results of these tests showed that the gain ratio 
of the two pairs of control ports did not remain constant as 
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^i~ure 4.1 

Typical Signer- Vipliiier Test Results 
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the control bias level or as the supply pressure was changed. 
However, for any given set of control bias pressures and a 
given supply pressure, the ratio of the gains of the two con- 
trol ports could repeatedly be obtained after signals were 
introduced and then removed. 

The results further showed tfhat f once the ratio of gains 
had been determined at a particular setting of bias levels and 

I 

supply pressure that the algebraic sum of two control signals 
could accurately be predicted. 

From the results of these tests, it tvas conoluded that 
this device could be used as a summing device with accurate 
and predictable results if the bias levels c^f the input sig- 
nals and the power ’'jet pressure were known. 



i 



I 



* 



62 



CHAPTER V 



DRIVING- BISTABLE AMPLIFIERS WITH PROPORTIONAL AMPLIFIERS 

1 

5.1 Introduction to, Input Control Impedances of Bistable 
Amplifiers ^- 2 

In Figure 5.1 a typical bistable fluid jet amplifier' is 
shown with the power jet attached to the left boundary wall. 

In this condition, the entrainment flow through the left con- 
trol port is greater than that through the right when both 
control ports are open to the atmosphere. If the left con- 
trol port is sealed from the atmosphere while the right port 
is left open and a pTdt of the left control pressure and 
flow is made, a plot similar to Figure 5.'2 will be obtained. 
This graph shows that as the left control pressure is in- 
creased above atmospheric pressure, a point is reached where 
the jet is supplied all the entrainment flow it can take. It 
then becomes separated from the left wall and switches to the 
right. When this occurs, the flow suddenly decreases cue to 
a different pressure in the interaction chamber and now pre- 
sented to the left control port. If the left control port 
pressure is then decreased, a point is reached where the jet 
will switch back to the left side. 

If the amplifier has perfectly symetrical control ports, 
the "stream-on-right” curve will represent the pressure-flow 
characteristics at the right control port when the stream is 
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r'irurc 5.1 

IstaMe fluid Jet --npli icr 




v iC'ire r.2 

Input Char no t eristics of a Ji stable -.rplifier 
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attached to the left wall. Prom the graph it is seen that the 
flow differential required for switching from left to right is 
and the pressure differential required is P A -? atm 
Similarly the flow differential to switch from right to left 
is 0,3-0,^ and the' required pressure differential is 
The two sets of flow and pressure differentials are not nec- 
cessarily equal. 

The effect of loading on the output legs of these de- 
vices is generally to decrease the flow and pressures required 
for switching. 

5.2 Problem Statement and Objective of Experiments 

Because of the nature of operation of proportional 
amplifiers, the useful signals from these devices are elevated 
about a bias pressure level greater than atmospheric . '• It was 
not clear what the effect of this bias level would be on the 
switching characteristics of a bistable amplifier being driven 
by a proportional amplifier. This information was necessary 

1 

before undertaking the design of the bistable reaction jets 
in order that the pressure differential required at the con- 
trols of the proportional amplifier, to effect switching of 
the reaction jets, might be minimized. 

The test procedure was to bias the control ports of a 
memory device' 11 at various pressure levels above atmospheric 
and then measure the pressures and flows required to effect 
switching at these elevated levels. 



* A memory device was selected for these tests as this 
is the type of device used for the first stage of the reaction 
jets discussed in the next chepter. 
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During the tests the output leg^-were loaded into an 
orifice whose area was approximately 0.75 times the receiver 
areas of the device. 

♦ 

♦ 

5.3 Test Hesults and Conclusions 

The memory device selected for these tests was a 

l 

Corning Standard Size .JLoad Insensitive device*. 

A graph of representative results is shown in Figure 5.3. 
It can be seen that the effect of the bias level is to radi- 
cally increase the flow and pressure differential required for 
switching. ' \ 

4 

■ During these tests it was found that the input charac- 
teristics of each control port were quite different. There- 
fore, to more accurately show the flow and pressure differ- 
entials required for switching, the pressure-flow relations 

♦ 

for both legs are plotted. 

Y.'hen a bias level is present at the control ports of 
a bistable amplifier, the power jet is being supplied more • 
flow than it can entrain. In this condition the bias flow 
is attempting to operate the device in a proportional manner. 
However, due to the close proximity of the bounding walls, 
the jet 'will still reduce the pressure slightly on one side 
and attach itself to that wall. This increased pressure 
differential required for switching can not be tolerated in 
this system as it represents a large hysteresis in the 



* Catalog No. FD 2212-2-1211 
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Effect of Biasing the Controls of a Bistable Bnplifier 
on the Switching Pressures and Blows 
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switching of the bistable reaction jets. This problem, how- 
ever, can be eliminated by making the bistable amplifier 
large enough so that the entrainment flow of the power jet 
is just equal to the bias level flow. 



$ 
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CHAPIER VI 



POWER AMPLIFIER AND BISTABLE REACTION JETS ANALYSIS AND DESIGN 



6.1 Determination of the Required Output Thrust 

So far little mention has been made of the dynamics of 
the system to be controlled. As mentioned in the introductory 
chapter, the system was to approximate an ideal inertia load. 

The dynamics of such a system are given by 

• • 

T = zQ (6.1) 

where T = the applied torque to the system 

J - polar moment of inertia 

0 = angular displacement 

• - differentiation with respect to time 

If this equation is integratea once with respect to 
time and the initial angular rate is assumed to be zero, one 
obtains 





( 6 . 2 ) 



Integrating equation (6.2) with respect to time and setting 
the initial angular displacement to zero gives 




T 



t 



2 



“ 2 “ 



(6.3) 



69 



Fro n the design considerations discussed in the in- 
troductory chapter, it was desired that the system respond 
fairly slowly and be able to receive large angular displace- 
ments. On the basis of these considerations, it was decided 
that the parameter T/J shoulc be such that if the system 
starts from rest under full acceleration', it must be travel- 
ing at an angular velocity equal to 1.1 radians per second 
at the end of 45 degrees of angular displacement. 

If equa.tion (6.3) is solved for t and substituted 
into equation (6.2) the result upon rearranging is 







Substituting the values of CO and Q into this equation gives 

T/J = 0.772 seq"^. From. equation (6.3) the time to travel 

« 

45°, with this T/J ratio, is 1.424 seconds. 

Having determined the T/J" ratio, an estimate of the 

t 

systems final polar moment of inertia was made by accurately 

determining the polar moment of inertia of the rate sensor and 

multiplying the result by a factor of three to account for the 

remaining equipment to be added to the system. The J so cal- 

2 

culuted was found to be 7.5 in-ozf-sec which results in a 
required torque of 5.79 in-ozf. If the lever arm is assumed 
to be 5.79 inches, then the thrust required is one ozf. 

> 

6.2 Introductory; Theory of staging Unvented Bistable 

, + A ov , 12 , 13, 14 

ivmnli Tiers ’ 

This system was to develop the required thrust by 
using the change of momentum in a reaction jet. Since the 
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power available from the output legs of the summer- amplifier 
was too small to develop the required thrust, this signal 
had to be amplified. 

Conpagnuolo^ has shown that this power amplification 
can be efficiently accomplished by staging bistable amplifiers 
in ascending order of size with each stage operating at the 
same supply pressure. He has pointed out that in order to 
make such a cascade, two design principles must be adhered to: 

1) For maximum pressure recovery, the splitter cis- 

tances must be kept small, preferably in the core region of 
the power jets which in general is approximately six nozzle 
widths down stream. : ' 

2) In order that the units be capable of being loaded, 
each stage must have a high degree of stability which implies 
moving the splitter down stream to over approximately twelve 
nozzle widths. 

At the present state of the art, a bistable unit can 
not possess both the characteristics of high pressure re- 
covery (high efficiency) and high stability. As the splitter 

is moved down stream, the pressure which can be recovered, 

* 

falls rapidly. If the splitter is moved up stream, then the 
ability to operate under load, ( stability ) falls off rapidly. 



* The core region of the power jet is that region in 
which the velocity of the jet has not been reduced from its 
value at the power nozzle due to viscous friction with the 
surrounding fluid. 
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Compagnuolo suggests that since the last stage is the 
most important for maximum efficiency that it he mace with the 
splitter located in the core region of the power- jet. This 
stage con then be held into the load by making the preceeding 
stages high stability devices. The result is that the overall 
cascade exhibits high efficiency and stability under load. 

In general the instantaneous flow gain of high stability 
devices is approximately ten. This means that it is possible 
to increase the size of the power nozzle of each stage by 

» 

roughly a factor of ten. It should then be possible to have 
flow gains in a three stage device of at least 1000 and, if 
the efficiency of the cascade is good, power gains of greater 
than 5000. 

6.3 Initial Experimental Design 

Figure 6.1 shows an output leg of a bistable amplifier 
turned through a 90 degree angle. If a control volume is 
drawn as shown and the momentum equation for steady state 
operation is written for the x-direction, the resulting 
thrust due to the exiting jet is found to be: 

T = PA e V e 2 (6.5) 

Assuming uniform velocity over the exit area, equation (6.5) 
can be rewritten as 

T = p (6.6) 

A e 

From this last equation, knowing the thrust desired and 
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figure 6.1 

Control Volume for Determining the Output Thrust 
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assuming a reasonable value for the exit area, one can deter- 
mine the flow required to develop the desired amount of thrust. 

If one further knows the differential flow required 
for switching the first stage, then the flow gain needed can 
be readily determined. 

It was shown in the preceecing chapter that it is de- 
sirable to have the entrainment flow of the first stage 
approximately equal to the output bias flow of the driving 

proportional amplifier. This entrainment flow depends on 

1 2 

the geometry of the amplifier and the supply pressure. It 
con be seen thet if sufficient experimental data were avail- 
able that it would be possible to 

1) determine the supply pressure required for a chosen 
first stage amplifier design from the entrainment flow versus 

I 

supply pressure characteristics and 

2) with the supply pressure chosen as above, deter- 
mine the differential flqw required for switching. 

It would then be a simple matter to determine the flow gain 
required end the number of amplification stages needed. 

Lacking any such quantitative data, a preliminary 
two stage bistable amplifier was constructed in order to 
obtain the following information: (listed also ore the rea- 

sons -this information was desired) 

1) Entrainment flows of a first stage design at vari- 
ous supply pressures to determine the supply pressure setting 
to match the output flows of the summer-amplifier. 

2) Stability under load of the unit to check the set- 
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ting of the splitter distances, wall set back distences, etc. 

3) Output pressure -flow relationships to determine 
the optimum loading for maximum torque output. 

4) 'The number of stages required to develop the 
desired torque. 

5) The thrust output for various supply pressures and. 

loads. 

The important dimensions of the amplifiers used 

y 

in this initial design are as follows:* 

• First £tage 

Power nozzle width * w s 0.020 in. 

Control nozzle width =• w 
wall setback ■* 3w 
Splitter distance a. 14w 
Aspect ratio = h/w =. 1. 5w 
Receiver width - 3w 

Second Stage 

* 

.Power nozzle width s w a 0.0625 in. 

Control nozzle width - w 
wall setcack = 1.4w 



* Many of the dimensions given here, such as wall set- 
back and receiver width, are based on carefully considered 
data given in references 12, 13, and 14. Time and space lim- 
itations do not allow a full discussion of the effects of 
varying these parameters. 
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Splitter distance a 6 w 
Aspect ratio = 1.535w 

Receiver width s 1.5w 

I 

The devices were made out of plexiglass using a pan- 
tograph engraving madhine. 

The results of tests performed on this initial design 
showed that at a supply pressure of two psi the entrainment 
flow in the attached side was approximately 0.04 SCFtt. This 
flow was approximately equal to the no load output flow of 
the summer-amplifier when operating at the supply pressure 
expected for final design setting. 

This device was quite ‘stable provided the output load 

4r * 

was relatively small. Typical output curves for the device 
are shown in Figures 6.2 and 6.3. In Figure 6.2 the device 
is not being driven by an external signal ana therefore can- 
not be heavily loaded. The overload area represents condi- 
tions under which the second stage power jet commences to oscil- 
late. In Figure 6.3. the device is being driven by another 
bistable amplifier. It is apparent that the device can now 
be more heavily loaded before the overload region is reached. 

In order to determine the value of the exit nozzle 
area for maximum thrust, a number of nozzle load curves were 
drawn on a plot similar to Figure 6.3. It was found that 
maximum thrust output of one nozzle corresponded closely to 
the point of maximum power transfer between the amplifier 
cascade and the nozzles. i>ince loading decreases the flow 
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out of one leg, part of the flow is diverted into the other 
leg where it develops thrust in the opposite direction. There- 
fore, considering the total thrust developed from the device, 
it was found that the point of operation for maximum total 
thrust ley slightly up the output curve toward the region of 
higher flows. 

The results also showed that in order for this two 
stage device to develop an ounce of thrust that the supply 
pressure would have to be over ten psi. Since this was in- 
compatible with the entrainment flov; desired, the final 
design would have to have at least three stages. 

6 .4 Final Three Stage Amplifier Design 



first two stages in the final design and alter th6m only 
slightly in order to accommodate the third stage. This 
further allowed the power nozzle area of the third stage to 
be calculated from the data obtained on the two stage device. 

The procedure used to determine the dimensions of the 
third stage are as follows: Assuming incompressible flow*, 

the flow from the third stage power nozzle is given by 



* This assumption is reasonably valid as the final 
device was to operate at two psi. 



In order to use the -data gained from the two stege 
amplifier to full advantage, it was decided to use these 




(6.7) 
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where Q,j3 = flow from third stage power nozzle 

a 53 = area of third stage power nozzle 

= discharge coefficient 



Similarly the second and first stage flows are given by: 



q j2 


: 0a 


- p s' 


(6.8) 


*ji 




p s 


(6.9) r 


The desired thrust is 


given by 
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(6.10) 



A j3 



where it is assumed that the load orifice areas are equal to 

the areas of the third stage power jet and the factor 0.8 is to 

* 

account for leakage flow out the opposite leg. 

From the two stage amplifier, numerical values are 
known for equations (6.8) and (6.9). Furthermore, since ? s 
is known and assuming a value for c^, equation (6.7) reduces 
to _ ' 

Q,j3 = (Constant) Aj3 (6.11) 

If the numerical values of equations (6.8J and (6.9) 
are substituted into equation (6.10) along 'with equationl 6.11) , 
a quadratic equation in A is obtained. Solving this 
quadratic equation gives the third stage nozzle area. Since 
all dimensions of the third stage are determined from the 

f 
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nozzle width, assuming an aspect ratio for the third stage 
will allow the nozzle width to b.e determined and hence spec- 
ify the other dimensions of* the third stage. 

Following the above procedure resulted in the following 
dimensions for the final design: 

First Stage 

Power nozzle width * w. - 0.020 in. 

Control nozzle width 1.5w^ 

. Wall setback = 3.0w^ 

Splitter distance = 18w^ 

Receiver width * 4w^ 

Aspect ratio » 1.5w^ 

2 

Nozzle area = 0.0006 in. 

Second Stage 

Power nozzle width 3 w 2 * 0.0625 in. 

Control nozzle width - w 2 
Wall setback » 2.5w 2 
Splitter distance **• 14w 2 
Receiver width * 3w 2 

Aspect ratio » 1.535w 2 

O 

Nozzle area = 0.006o5 in. 

Third Stage 

Power nozzle width =» w^ » 0.125 in. 

Control nozzle width - Wj 
’Wall setback 3 2w^ 
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Splitter distance 
Receiver width 
Aspect ratio 
Nozzle area 



6w 3 
1.4v/ 3 
2.08w 3 
0.0325 in. 2 



A photograph of this device is shown in figure D.4. 

The circular chambers are vortex impedance matchers 
modeled after those used by Compagnuolo^^. 

In Figure 6.4 is shown the output thrust of the ampli- 
fier cascade versus supply pressure and input flow rate. The 
output flow of this final design oscillates due to the load- 
ing conditions unless the cascade is being driven by an ex- 
ternal signal, 
l 

Figure 6.5 shows the input impedance characteristics 
for the right control port when the first stage power stream 
is attached to the right wall. The supply pressure is 1.96 
psi. 



6 . 5 Reaction Jet Cascade Dynamics and Transfer Function 

The reaction jet cascade is characterized by two 
phenomena which are important when the device is placed in the 
control loop, of the final systeiij: 

1) A hysteresis loop is ^present in the switching of 
the thrust from one side to the other 

2) There is a finite time delay between the time a 
switching signal is received at the control ports until the 
thrust developed is switched from one side to the other. 

4 

The hysteresis loop is due to the characteristics of 
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figure 6.4 

Reaction Jet Thrust Versus Supply Pressure 
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} Reaction et Input Characteristics 
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bistable devices requiring a finite flow and pressure dif- 

/ 

ferential to exist across the control ports to effect switch- 
ing. The half width of the hysteresis band was determined 
from actual measurements with the summer-amplifier set at its 
proper supply pressure and driving the reaction jet cascade. 
The half width was found to be 0.0685 psi. The time delay 
involves two major phenomena: 

1) The time it takes the fluid streams to traverse 
physical distances (transport time) and 

2) The time it takes to charge the interconnecting 
volumes to the required switching levels. 

Actual measurements of this time delay were not made. 
From theoretical calculations in Appendix B the time delay 
constant was found to be 0.022 seconds. 

The transfer function of the reaction jet cascade is 
thus given by 

* 7 * 3 

± T - e 3 with hysteresis oft 0.06e5 

£P psi 

( 6 . 12 ) 



85 



I 



CHAPTER VII 



CONTROL SYSTEM SYNTHESIS , PEREC PJiANCE AND ANALYSIS 



7.1 oyster. Synthesis 



A major part of the problem of synthesizing the system 
components has been accomplished by cereful consideration of 
the role each of these components was to play in the final, 
system. However, in order to finally bring these components 
together into a working system it is necessary to know hov< 
each component will behave when actively loaded into the 
component it is to drive in the final design. 

Of critical importance is the setting of the power 
supply pressure to the summer-amplifier. This pressure must 
be high enough to prevent the reaction jet cascade from oscil- 
lating and, at the same time, allow for maximum amplification 
of the rate and position control signals. It can reacily be 
seen that had the reaction jet cascade been designed with no 
consideration for this supply pressure setting, that the mag- 
nitude of the control signal required for switching the jets 
might have been so large. as to preclude the operation of the 
system. 

The rate sensor had been tested with its outputs loaded 
into passive orifices approximating the control orifices of the 
summer-amplifier. It was thus necessary to recalibrate the 
rate sensor with its outputs* loaded into the summer-amplifier. 



86 



For accurate results the summer-amplifier must be set at its 

proper supply pressure and, ,in aduition, must have the second 

set of controls activated . 1 This is necessary as the input 

impedance of the control nozzles is slightly affected by the 

supply pressure and bias level at the other controls. The 

steady state results of these calibration tests on the rate 

sensor are shown in Figure 7.1. For these tests, the rate 

« 

sensor output was loaded into the summer-amplifier with the 
power supply pressure set at its operating level. The second 
set of controls were not activated. The results shown in this 
figure were not obtained from an accurate rate table and are 
considered reasonably accurate only in the range ^ 0.5 radians 
per second. 

In order to set the maximum pressure signal from the 
position sensor*, it was necessary to know: 

1) The gain ratio of the twc summer-amplifier control 
ports and 

2) The final system T/J ratio. 

As was pointed out in the discussion of the signal comparator, 
the gain ratio depends on the summer-amplifier power supply 
pressure setting end control bias levels. Thus, to obtain the 
correct gain ratio, a new set of curves similar to those 
shown in Figure 4.1 were constructed for the proper supply 
pressure setting and control bias levels. 



* The maximum signal occurs at ±90 degrees of angular 
rotation. 
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Pifure 7.1 

Steady State Angular Rate Sensor Output 
•vhile Loaded into the Sunner— rplifier ' 

l 
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As shown in Appendix C.lJ the T/J ratio defines the 
switching curve equation. The criteria used to sot the maximum 
pressure signal from the position sensor was as follows: If 

* i 

\ 

the system starts with an initial error of Q z 90 degrees and 

# 

Q - 0 radians per second, the maximum rate signal at the 

k 

point of theoretical optimum switching will occur at Q - 

•i 

45 degrees. To effect switching of the jets, the position 
signal should then equal the rate signal ot this point. The 

0 

T/J ratio defines the value of this peak^ . From Figure 7.1, 

# 

the pressure signal corresponding to this*value of Q is deter- 
mined. This then is the value that the position signal must 
have at 45 degrees. Knowing the shape of the real switching 
curve and the value of the gain ratio of the summer-amplifier, 
the value of the maximum position signal can readily be deter- 
mined. * 

In the final system, a pair of adjustable orifices 
were placed between the output of the position sensor and 
the input to the summer- amplifier to make fine adjustments of 
the switching curve. 

Figure 7.2 shows a block diagram pf the final control 



system composed of the 


various 


transfer functions from the 


previous chapters. A 


summary 


of the numerical values for 


the parameters is given below: 




s. 


1.19 


. psi 


X 


rad/sec. 


*K 0 s 


1.14 


psi 


2 


rad 



*These gains have included the gains of the summer- 
amplifier. The summer amplifier multiplied the rate signal 
by six and the position signal by 4.5 
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Figure 7.2 



Block Diagram of Complete Control Syste: 



0.055 sec. 



7*2 = 0.04 sec. 

T 3 = 0.022 sec. 

T - 7.72 in-ozf 

hysteresis* —0.0666 psi 

J - 9.60 in-ozf-sec 2 

Figure 7.3 shows thei same system in circuit disgram term*. 

A photograph of the complete control system is shown in 
Figure 7.4. 

7 . 2 System Performance 

To obtain performance data on the system, a rotary 
potentiometer was connected to the control system and the 
output signal displayed on an X-Y recorder. The results are 
shown in Figures 7.5 thru 7.9. 

Figure 7.5 shows the response of the system to a num- 
ber of position errors. The reaction jets switched in a near 
ly perfect manner when responding to counter clockwise errors 
The response to clockwise errors, however, resulted in a 
slight amount of "after end point chattering". The magni- 
tude of the limit is seen to be erratic, averaging roughly 
8.5 degrees at approximately 0.5 cps. This erratic response 
is due to noise from the proportional amplifiers in the 
system. Considering the T/J ratio of this system, the 



* The symbols used in this figure follow those reconv 
mended by Boo the l- 5 . 
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figure 7.3 

Systen Circuit Diagram 



I 
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Figure 7.4 

Complete Control System and Demonstration Model 
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response time for the errors shown was found to be nearly 
optimal. 

Figure 7.6 shows the system operation with a slightly 
lower summer-amplifier supply pressure setting. By lowering 
the supply pressure, the magnitude of the hysteresis oar.c was 
decreased, resulting in the marked decrease in the magnitude 
of the limit cycle. Unfortunately, the noise present in the 
system could then switch the reaction jets more readily re- 
sulting in a slightly less optimal time response to large 
errors because of random switching. 

Figure 7.7 shows the system performance with the rate 
sensor disconnected from the control loop. As pointed out 
more clearly in Appendix C.2, the control system should, 
theoretically, be completely unstable. Because of the finite 
camping in the system, in practice a new limit cycle was 
reached at a much higher emplitide and lower frequency. 

Figures 7.8 and 7.9 show an interesting mode of oper- 
ation, possible with this system, ns pointed out in the pre- 
ceding chapter, the reaction jet cascade is unstable unless 
driven, due to the heavy loading on the: output legs of the 
last stage. This back pressure is so high that it causes the 
first and second stage amplifiers also to oscillate. Intro- 
ducing small control flows into one side of the first stage 
cause the pulsing to occur more often on one side than on the 
other in direct proportion to the magnitude of the input 
signal. Thus the system effectively operates in a pulse 
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width modulation mode * . 

In figures 7.S and 7.9 this pulse width modulation node 
was obtained by lowering the summer -amplifier supply pressure. 
The frequency of oscillation is too high (approximately 
100 cps) for the system to respond to the individual torque 
pulses . 

r 

Figure 7.8 shows the response ,to various disturbances 

« 

while operating in this pulst width modulation mode. Of 

particular note is the disappearance of random fluctuations 

due to noise. A certain amount of stiction is noted caused 

• » 

! 

by the potentiometer attached to the system. 

Figure 7.9 shows the system response with the rate 

; * 

sensor disconnected from the control loop. If the load were 

•i 

completely free of camping, the response, of course, would 

not camp out. A comparison of this figure with Figure 7.7 

shows the functioning of the vortex rate sensor and ttye part 

it plays in the stabilization of the system. 

* 

7 . 3 System .-.r.e I'-s 1 s 

This control system falls into the category of non- 
linear control systems because of the cang-bang method of 
controlling the load. The most convenient way of analyzing 
this system, in order to predict the magnitude and frequency 
of the limit cycle and the stability of the system, is by 
using the cescribing-f unction method. 

Gibson^ points out that this approximation method is 



*For a more complete explanation of this mode of oper 
at ion see reference 17. 
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quite accurate in systems containing sufficient lev; pass filter- 
ing. The nature of an inertial lose amply fulfills this re- 
quirement as it refuses to respond (filters) to high frequency 
signals . 

This method essentially consists of assuming that yhe 
input to the non-linear element is sinusoidal. The resulting 
output is then written in the form of a 1‘curier series. Then 
it is assumed t^at only the fundamental component of this 
series is important in relation to the response. The coeficient 
of this fundamental frequency, which is dependent on the mag- 
nitude of input signal but independent of the frequency, is 
used as an ecuivalent gain (K_„) of the non-linear element. 

The rest of the system being linear is frequency dependent 
but not dependent on the magnitude of the input signal. By 
plotting the equivalent, gain function of the non-linear 
element and the frequency response of the linear elements on 
a polar plot, a point of intersection may be found which re- 
presents a possible point of operation. If the point is 
found to be a staole point, then the frequency and magnitude 
of oscillation can readily be determined. 

This method is applied to the analysis of this system 
in Appendix (J.2 The results closely predict both the magni- 
tude and frequency of the limit cycle shown in Figure 7.5*. 



*The data on which the calculations were based was 
obtained from the system just after Figure 7.5 was recorded. 
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APP2NL1R A 



RATS SSNSCR CALCULATIONS 
A . 1 Rate Sensor Design Calculations 







-« 


Given 


*0 = 


3 in. 




k^nax - 


70°/sec =• 1.221 rad/sec 


A ? *s 


CO = 

max 


0.1 lbf/in 2 


A(5> w 


^ max 1 


3° = 0.0523 rad. 




P • 


2.34L10" 3 — -^'-4 - 




y r 


1.623L10” 4 ft 2 /sec 



From equation (2. IS) 

- “\ / 144 A~? 
’ V 4C 



u c 



4 P A 0 



-v / (144) (0.1) 

V / 



4(0.00243) (0.0523) 



U s - 171.2 ft/sec 



From equation (2.21) 



r\ — 

^ * 



^r s - 4cr s 



U 



[^(0.1562) 2 - 4(0.0212) (0. 1562)] 171.; 



144 



Q, * 0.151 ft 3 /sec 9.05 ft 3 /min. 
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From equation (2.22) 



Re = 

y P 2 



He s (0.151) (0.0251) 

(1.623a15 4 ) (0.0625) 
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From reference 3, page 327, vje find that, this value of He is 
off the scale but E is leveling off to roughly 0.66. k 
"factor of safety" was introduced here and E was taken to 
be 0.5. 

Rearranging equation (2.20) gives 



Ro 2 k) 2 

r P° 1 A$ u s 



r 



po 



(0.0625) (1.221) (0.5) 
(0.0523) (171.2) 



4.27 a tcf^ft 



-2 . 

5.12 a 10 , me. 



po 



Rearranging equation ( 2 . 23 ) gives 



a 



a 




a 



0.0475 in 



k .2 Calculation of First otuge oupoly Pressure 



From Figure 2.8, Q, 



c max 



i 0.107 in 3 /sec and take Q, c 0 



Then 




0.1C7 
3. 75X10”"* 



285 in/sec 



10 ^- 



fror. equation (2.24), V D = 



v ; c V 

w^ tang) 






where w p/ 2 



0 .005 
0.1 



= .05 



w c = 0.015 in. 



Wp = C.01C in. 



Therefore, V p 2 2.43 jC 10 6 in 2 /sec 2 



From equation (2.28) 



‘P V P = V C cp ^P 




or 



V 2 

P D 

P 



c dp 2 



P/2 



Fron. Corning design data sheets, C^ p =. 0.70, 



Thus P p *0.28 psi - 7.7 inches of v; ater. 



i 
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APPENDIX B 

* 

Gr.LGULitTIOII Or THE RE^CTiw* JET Till- GELAY CONSTANT 



The time required for the power jet to traverse the dis- 
tance between the nozzle and the receiver ports is given by 

t = IT— 

. avg 

v;here 1 - splitter distance 

V avo . - the average velocity of a particle 
while traveling between the power 
nozzle exit and the receiver. 

Reference 19 gives the dynamic pressure profile of a sub- 
merged jet at various distances from the power nozzle. If 
it is assumed that the velocity decreases linearly from noz- 
zle exit to the receiver, then V is given by 



nozzle exit velocity 

V s - jet velocity at the entrance to 
the splitter 

i 

The value of V 3 can be determined from the grephs given in 
reference 19 if the splitter distance is known. For the first 
stage : 

splitter distance * 13w s 0.36 in * 1 






avg 



h-ll s 



where V 



j 



■V 



T ^ 



t 
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iaLJJU2£i 

l i. 129 a 10“') 



5.9 X 10 3 
in/sec 



Ell Via2 = 0.45 
P/2 V 0 - 2 



V I3 c v s = 3 » 96 A 10 3 in/sec 



avg 



v ,i y s = 4.93 a 10 3 in/sec 



.t- 



_ 0.35 



V V 



avg 



4.93 .v 10' 



= 73 a 10" 6 



secs. 



Similarly for the second and third stages: 



t’ 2 = 197 & 10~ 6 • secs. 

tr. - 137 a 10" 6 secs, 

o i 

Che time required for the- flow in the control nozzle of 
an intermediate stage to build to that flov; level necessary 
for switching can be approximated from the time constant as- 
sociated 'with the volume connecting the two stages and the 
resistance of the control port. The time constant is given 
by 

*T - RC 

where R = impedance of the control sort lbf -sac/in 5 

« 

C : capacitance due to the interconnecting 

volumes in 5 /lbf 

The value of R was determined by linearizing the input impedance 
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curves of each, nozzle about an assumed operating point. The 
value of C is given, .for adiabatic conditions, by 



C = 



where v - 

P- 

K = 
& = 



v 

p icftr 

connecting volume (in' 5 ) 
mass density ( Ibf-sec^/in^) 
adiabatic constant = 1.4 

gas constant = 2.48 JL 10^ ( in^'/sec*’-'*^) 
temperature (°k) 



For lines connecting the first and second stage: 

1 

R = 0.287 l'of-sec/in^ 

I 

The value of v was found to equal 0.279 in 3 . Therefore 

I 

c _ v _ (0.279) 

P&R T " ( 1.129X10?) (1.4) (2.48X10 b ) (550) 

C = 0.013 in 3 /lbf 



Eence 



T - .00372 secs. 



The relationship between the flow into the connecting volume 
(assumed constant") and the flow through the control nozzle 
(variable) is given by 



^in 7" S4*l 

The value of is known from the nozzle geometry of the 
first stage. ^ o i- 3 zero before switching occurs. Af t er 
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switching *; 0 increases until it reaches the value required 
for switching. The value required for switching was taken 
to be nearly equal to «iin* Therefore, the tine required for 
the flow to build to the switching level is approximately 
four tine constants, or t 0.0140 seconds. In a similar 
manner the buildup tine between the second and third stage 
was found to equal 0.006 seconds. The final time delay is 
due to the tine required for the output flow to traverse 
the length of the output leg. By an elementary calculation 
this was found to equal 636 L 10“° seconds. u 

The total tine delay is equal to the sun of all the 
tine delays discussed ebove, 0.022 seconds. 
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APPENDIX 



SYSTEM ANALYSIS 



C.l Determining the Equation of rhe Switching Curve 



Consider^ a bang-tang control system as shown in figure 
C.l. The equation of motion describing the operation of this 
system is given by 



& 



Sign Q 



T_ 

J 



(c.l) 



if 



e 



is rewritten as 



Q- dCJ T dO) 60 . ^ d CO 
dt d Q dt d^ 

and substituted for Q in equation C.l one obtains 

sign^> 



T_ 

J 



CO dCO 
! ^0 

Equation C.2 can now be integrated resulting in 



( C . 2 ) 



(o^ . _ _£_0 -rCq for Q (C.3a) 

2 ’ J 

_ T C\ Cg for ( + ) Q (tC.3b) 

2 J 

These curves represent two families of parabolas in the phase 
plane (CJ vs 0 ) where Cq and C 2 determine the values of the 
initial conditions. If these curves are plotted in the phase 
plane, it will be found that the curves corresponding to 
C]_ =• C£ *0 pass through the origin. They represent trajectories 
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figure C.l 

Idealized Contactor Servomechanism Controlling 

an Tnertia Load 




Eigure C.2 

Generalized Autonomous Control System Containin' 7 
a Non 1 inear Element 
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such that, if the system. is switched at the tine one or the 
other curve is intersected, a tine optimul reduction of errors 
introduced into the systen will result. 

Thus the switching line equations are gjLven by 



CO 



9 



2 




second quadrant 
switching curve 




fourth quadrant 
switching curve 



G . 2 Limit Cycle Analysis 



The bloch diagram shown in Figure 7.2 can be reduced 
to the form shown in Figure G.2 if the following conditions 
exist: 

1} The reference signal is zero ' 

2) The magnitude of the position signal stays 
within the linear region of -the switching curvb.* 

In the figure, K represents the nonlinear element in the 
system and G is the combination of ail the linear elements. 

Gibson 1 shows that a sufficient condition for the 
existence of sustained oscillation is 

G( joJ) 

Keq 

where Keq =. the equivalent gain of the non- 

linear element, which is a func- 
tion of the magnitude of the input 
signal to the nonlinear element. 



* This condition is satisfied os shown by the results 
in Figure 7.5. 
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Eeq for a contactor with hysteresis is given by Gibson as 



Keq = g (2) + j b(E) 



where 



and 



2 = the magnitude 'of the input to the 

contactor represented by Aly.-*- A 

U U 



g(E) -.-£L 

' irs 



V 1 - ( '2-> 



b(E)r-4T_ _R_ 

'ITS E 



whe re 



T 



the half v;idth of hysteresis band 
(psi ) 

the magnitude of the torque output 
( in-ozf ) 



For this system, 



0.0686 psi 



7.72 in-ozf 



By ^ssuning values for S*A^+-A^ , various values of - l/Keq 
may be obtained. For example 

assume E = 3 inches of water 0.1082 psi. 



g 5 



g = 




1 - 



4 (7.72) 
-nK 0.1032} 



* 

V I -/^ 0«06G6 

\0.1032 



* These values are normalized Fourier coeficients. 
The term b(E) is needed to take care of phase distortion of 
the fundamental harmonic. 
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s ? 



70 



' iP-ozf 

psi 



h - 4Te - _ 4(7.72) (0.0586 _ _ = 7 , in-ozf 
'TTES ' (0.1082) ^ ~ * psi 

Keq - 89.5 > r - 39.4° 





, 1 . 


0 01 - 2 


psi > 


39.4° 




Xeq 


V t V + 


in-ozf ‘ 


The linear 


transfer 


function 0(s) is given by 












t 


G(S)- 


^Se” 


Tjs^ :< s 


e-T 2 3 






_ 


r 2 s*i j 










JS 2 




Assume that 


the time 


delay 


of the rate 


sensor 



(C.4) 



zero in order to see whet effect rate feedback has on the 
system. For this case, the linear transfer function becomes 

-r 3 s 



Cr( 3 ) g ^2 






(C.5) 



S 2 (T 2 S*l) 



3y substituting actual numerical values gives 



I 

0(6) z (0.1188) 



S 2 + 2(2.55) 3+1 
(4.S9) 2 4.89 

2 

8 (0.043+1) 



-0.0228 

e 



( C.6 ) 



To show the effect of rate feedback on the stability of this 
system, disconnect the rate f.eec.back (X-^-o) . For this case 
equation (0.6) becomes 



- 0.0225 

0(6) - 0.1188 

S 2 (0.048 «#■!) 



11 k 



By substituting s-.jco, the values for the magnitude and phase 
of G(jdj) as CO varies from zero to infinity are obtained. If 
-his frequency response date is plotted on a polar 1 plot along 

with the magnitude and phase of , a plot similar to that 

shown in figure C.3 will be obtained. It can be seen that no 
intersection of the two plots occur end hence the system is 
unstable. 



Substituting s-jco into equation (C.6) and plotting this 
function on a polar plot (Figure 0.4) shows the stabilizing 
effects of rate feedback. It can be seen that a point of 
intersection now occurs and a stable limit cycle results. 

The effect of rate feedback is to introduce phase lead into 
the system which "pulls" the G(jdO) curve away from the -180 
degrees phase shift line a sufficient amount to allow this 
intersection to occur. To investigate what effects the time 
delay in the rate feedback will have on the system, Equation 
(C.4) may be written in the following form 



G(S) 5 K 2 



K n T. 



1 2 3 2 



-T ■, 3 



oe 



+ 1 



*2 



“2 



3^ (TpO+1) 



If the rate sensor time delay tern is factored out of the 
quadratic term and combined with the reaction jet time delay, 
the result is 



G(3) ; ^p - _ ^ B 2 -*- 

~f-[ X 2 



/ 

° ( r 2 s + i) 




e “l‘T 1 vr 3 j3 



( C . 7 ) 
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Figure C.3 




Figure C.4 

Polar Plot of G(jco) and - -1- .Vith Rate Feedback 

^eq 
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effect to be noted is the increased transport tine delay given 
by the tern Z ^ s . The other effect to bo noted occurs 



in the quadratic tern'. If the quadratic terns given in these 
two equations ere piotteo over the frequency range o 4.60400, 



by a very slight a noun t around the break frequency. Thus, the 
nsjor effect of the rate sensor tine delay is to increase 
the transport tine delay of the over all systen. 

Figure C.. 5 shows a Bode plot of equations (C.5) and 
(C.7) where the effect of the rate sensor tine delay can be. 
seen nore clearly. 



The effect of the rate sensor tine delay is to increase the 
nagnitude of the limit cycle and decrease the frequency. 

Fron this figure, the frequency of the Unit cycle for the 
systen is predicted to be 9.5 radians per second or 1.51 cps. 



wares through the linear transfer function to get an expres 




Figure C.6 shows polar plots of G(-'c 0) and 



nee 





This procedure gives 




where at the point of intersection 



o 



0.1G32 psi 
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nprnxxiu-'V 



I 




: \r, ire 0.5 

r o8e Plots of Linear Trans/'er ■'’’notion GCjco) 
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Fhase ^tift (degree*) 




i 




Figure C.6 

Polar Plot of C7(joj) and - 

^oq 
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CO 0 = 9.5 radians per second. 

Making these substitutions and carrying out the indicated 
calculation gives 



6 



max 



5.68 degrees 
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APPENDIX D 

PHOTOGRAPHS OF COMPONENTS 



T 



Figure D.1 

Vortex Rate Sensor Showing Details of Signal Amplification 





Figure D.2 Sumner -Amplifier 
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Figure D.3 

Angular Position Sensor 




Figure D. ! f 

Pover Anplil'ier Cascade and Reaction Jets 
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